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Abstract 



The spectrum of charged particles hopping on a kagome lattice in a uni- 
form transverse magnetic field shows an unusual set of Landau levels at 
low field. They are unusual in two respects: the lowest Landau levels are 
paramagnetic so their energies decrease linearly with increasing field magni- 
tude, and the spacings between the levels are not equal. These features are 
shown to follow from the degeneracy of the energy bands in zero magnetic 
field. We give a general discussion of Landau levels in the case of two degener- 
ate bands, and show how the kagome lattice tight-binding model includes one 
special case of this more general problem. We also discuss the consequences 
of this for the behavior of the critical temperature of a kagome grid super- 
conducting wire network, which is the experimental system that originally 
motivated this work. 
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The behavior of electrons in a two-dimensional periodic structure immersed in an exter- 
nal magnetic field has been a subject of special interest over the past two decades ever since 
Hofstadter's graphical solution for square lattice 1 . The competition between the periodic 
potential imposed by the lattice structure and the length scale provided by the magnetic 
field leads to frustration in the system and to intricate, detailed, and self-similar features 
for the subband spectrum of single-particle eigenenergies. The magnetic-field-induced frus- 
tration in two-dimensional structures has also been a subject of active experimental study 
using superconducting wire networks 2 ' 3 , as the linearized Ginzburg-Landau equation for the 
superconductor and the Schrodinger equation for the tight-binding electron wavefunction 
in these lattices share the same form 4 . The investigation of a variety of structures, such 
as the triangular lattice 5 , the honeycomb lattice 6 , and recently the T% or dice lattice 7 ' 8 , has 
revealed the richness and the beauty of physics in these systems. 

In this paper, we present our study of the Landau subband spectrum of the kagome 
lattice, whose distinguishing degenerate ground state has attracted extensive attention over 
the last decade. The kagome lattice is a two-dimensional lattice of corner-sharing triangles 
(Fig. 1(a)). The nearest-neighbor-coupled Heisenberg antiferromagnet system on a kagome 
lattice has been shown to have an infinite number of classical ground states, which give rise 
to a zero temperature extensive entropy 9,10 . However this degeneracy may be partially lifted 
by quantum or thermal fluctuations 11-13 . This so-called order from disorder mechanism 14 
selects the states with the largest fluctuations (thus the highest entropy or the lowest zero- 
point energy), which often are more ordered states. 

Experimentally, kagome systems have been previously studied in the context of adsorbed 
3 He on graphite at millikelvin temperature 9,15,16 , and the layered oxide SrCrs- x Ga4 +x Oi9 10,17 . 
Our recent experimental studies 18 of nearly perfect kagome structures using superconducting 
wire networks show interesting and complex phase boundaries between the normal and 
superconducting states, which agree well with Y.L. Lin and F. Nori's mean-field calculation 
based on quantum interference 19 . This work shows that a frustrated ground state exists at 
a magnetic field of one half of a flux quantum per triangular plaquette. 

The lowest energy states of a charged particle in a nondegenerate band become Landau 
levels when a small magnetic field B is added. The energies of these Landau levels increase 
linearly in \B\; the currents in these states are dmmagnetic. The consequence of this for 
the critical temperature T c of a periodic superconducting grid is a local maximum in T c at 
certain rational values of the magnetic flux per unit cell (including, of course, zero magnetic 
field), with T c decreasing linearly in the magnitude of the change of the field away from 
these special values, as has been observed in experiments on various types of wire grids and 
Josephson-j unction arrays. For the frustrated kagome grid, on the other hand, there are two 
bands of lowest energy states, and the currents in the resulting lowest Landau level of the 
tight-binding model are instead paramagnetic and the state's energy decreases linearly in 
\B\. This results in a minimum in T c vs. field at the frustrating value of one-half of a flux 
quantum per triangular plaquette, as is seen experimentally in the kagome superconducting 
grids at this frustrating value of the applied field. 

Thus from the specific kagome lattice tight-binding model, we are led to consider the 
general case of Landau levels in a small magnetic field with two degenerate bands. In this 
paper we examine this more general problem first, and then show how the kagome lattice 
spectrum reduces to a special case, which can be readily solved analytically. We have also 
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obtained extensive numerical results for the kagome tight-binding spectrum, which agree well 
with our analytic results. We conclude with a discussion of the kagome grid superconducting 
wire networks that motivated this work. 

I. TWO-BAND LANDAU LEVELS 

Here we will explore a generalization of the usual Landau levels for a charged particle 
moving in two dimensions in a uniform magnetic field. The new case we consider is where 
the particle in zero field has two bands that are degenerate at zero momentum. This is the 
case in the kagome lattice tight-binding model. 

What are the "rules" we will use to generalize the usual one-band case? The particle 
is moving in two dimensions. We look in the continuum limit of low momentum p, so only 
look at the energy to order p 2 . At this order in p, the spectrum of the particle's states in 
zero magnetic field (B = 0) is assumed to be isotropic in momentum space, reflecting an 
underlying system that has at least three-fold rotational symmetry. The system is time- 
reversal invariant for B = and the spectrum depends only on \B\, not on the sign of B. 
(Because the motion is two-dimensional, only the normal component of the magnetic field 
enters.) 

For a single band, the only Hamiltonian that satisfies all these constraints is the usual 
H = n 2 /(2m*), where 

n = P -^A (1.1) 

c 

is the gauge-invariant momentum. In a small uniform magnetic field, this gives the usual 
spectrum of equally-spaced Landau levels E n = ftw c (n + |), with oj c = \qB\/(m*c). The 
currents are diamagnetic for these standard Landau levels. 

If instead we have two bands that are degenerate at p = for B = 0, then in addition 
to the momentum, the single-particle states are labelled by a band index, which may be 
treated as a spin-1/2. The general Hamiltonian in this case consists of Hermitian 2x2 
(Pauli) matrices operating on the band index and Hermitian quadratic combinations of 
and Uy. In addition to the zero-angular momentum combination n 2 that is allowed in the 
one-band case, there are two Hermitian operators quadratic in n with angular momentum 
two: n^n^ + n^n^. and n 2 , — n 2 , as well as the magnetic field itself, which is proportional 
to the operator combination ^H^n^ — n^Hj). 

We have a number of constraints that restrict our choice of Hamiltonian; let us first look 
at the simpler case of B = 0. The general Hermitian Hamiltonian that is quadratic in p is 

tf = f>£,(p), (1.2) 

j=0 

where the Oi are the Pauli matrices, with <7 the identity matrix, and the Ei are real quadratic 
functions of p. This is readily diagonalized, with the resulting eigenenergies 

E±(p) = E ±^E 2 + E 2 2 + Ei . (1.3) 
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Our restriction that the B = spectrum is rotationally invariant then dictates that E (p) = 
d p 2 and 

E\ + El + El = dip 4 , (1.4) 

with <i , c?2 real constants. Another way of writing the Hamiltonian which explicitly breaks 
it into its parts with orbital angular momentum quantum numbers m = 0, ±2 is 

H = d a p 2 + z • op 2 + (x - iy) ■ a(p x + ip y ) 2 + (x + iy) ■ a(p x - ip y f , (1.5) 

where x, y and z are real 3- vectors in the space spanned by the 3 Pauli matrices, o"i, 02, 03. 
Eq.(1.4) requires that these three vectors form a mutually orthogonal triad, with x and y 
being of equal magnitude. Thus we may choose the basis for the band index so that this 
triad x, y and z are parallel to the 1, 2 and 3 axes, respectively. This leaves a three real 
parameter family of B = two-band Hamiltonians: 

H = d a p 2 + d 3 a 3 p 2 + di(cr_(p x + ip y f + cr + (p x - ip y ) 2 ) , (1.6) 

where a± = (01 ± icr 2 )/2 are the usual raising and lowering operators. The spectrum is 
rotationally invariant, as we required: 



E ± (p)=p 2 (d ±^d 2 + d 2 ) . (1.7) 

An important observation about this Hamiltonian Eq.(1.6) is that it conserves a total angular 
momentum, provided that we define the particles in the bands to have an "internal" angular 
momentum whose quantum number m ; = ±1 is the eigenvalue of <t 3 in this basis. [We will 
see below that such an "internal" angular momentum rrii — ±1 within each unit cell does 
indeed exist for the low-energy, low-momentum states of the kagome lattice tight-binding 
model.] The operators a± then raise/lower the internal angular momentum by 2 quanta, 
while (p x ± ip y ) 2 raise/lower the orbital angular momentum by the same amount. Note 
that the operators (p x ± ip y ) are not simply the usual orbital angular momentum raising 
and lowering operators L±; in addition to changing the orbital angular moentum they also 
change the radial wavef unctions. If neither band has negative energy states at B = 0, so 
the ground state is at E — 0, this corresponds to d > ^d 2 + d\ > 0. 

Now when we add the uniform magnetic field, the momentum operator p becomes the 
gauge-invariant combination II = p — -A. For B 7^ the two components of IT are noncom- 
muting variables with a c-number commutator proportional to the magnetic field: 

[U x ,U y }=t^. (1.8) 

Thus we can make harmonic-oscillator-type raising and lowering operators a* and a such 
that 



n, = ffV + a ) (1.9) 



and 



n, =*sign(gW^|iV -a) . (1.10) 
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Since the magnetic field itself is quadratic in II (Eq.(1.8)), a term linear in B may also be 
added to the Hamiltonian. This term may not contain the identity Pauli matrix, or it will 
generate a spectrum that is not invariant under inverting B, but it is otherwise unrestricted. 
Thus we end up with the general two-band Hamiltonian satisfying all our "rules" : 



H = (d a + d 3 a 3 )U 2 + di(<7_( n z + iU y ) 2 + a + (Tl x - m y ) 2 ) + Bb ■ a , 

where b is a real 3-vector. 

At this point let us note that the current density operator for this Hamiltonian 

contains, in additional to conventional-looking terms, a term proportional to d\ that is off- 
diagonal in the band index and linear in the gauge-invariant momentum n. It is this term 
that allows Landau levels involving both bands to exhibit paramagnetic currents. These 
currents arise from an "interference" between the amplitude in one band and the gradient 
of the amplitude in the other. 

The transformation between the n operators and the raising and lowering operators a 
and a) is singular at B = 0, so it is simplest to just consider one sign of qB at a time. In 
fact it is for qB < that the raising and lowering operators a and a* also raise and lower the 
orbital angular momentum, so let's do that case. Since we are motivated by superconductors 
we will take B > and q = —e* = — 2e < 0. This yields 

hp* r >ftp*d^ 

H = B[—(d a + d 3 a 3 )(2a^a + 1) + t^l(a^f + a + a 2 ) + b • a] . 

(1.13) 

When the 1 or 2 components of b are nonzero, this introduces terms that flip the band index 
without changing the orbital angular momentum, so do not conserve total angular momen- 
tum. We have not found a closed-form solution for this general case, but for the kagome 
lattice problem we are considering, these terms are not present. In this case, b\ = b 2 = 0, the 
total angular momentum is conserved. The orbital angular momentum quantum number 
m is also the eigenvalue of the number operator a^a and is restricted to be nonnegative for 
the case qB < that we are now considering. Thus for total angular momentum m of or 
— 1, there is just one state, with internal angular momentum quantum number = —1 and 
m = m + l. For uii + m = m > 0, on the other hand, there are two states, with to* = ±1. 
If bi — &2 — so that m is conserved, the Hamiltonian then reduces to a simple 2x2 matrix. 
The resulting Landau levels are in general not equally spaced in energy, and even when there 
are no negative energy states at B — 0, some of the Landau levels may be at negative energy 
for B 7^ 0, and thus paramagnetic. Now will look at the kagome tight-binding model, and 
show that in the low-momentum limit it realizes this two-band Hamiltonian in the simple 
case of do = d\ > 0, d 3 = 0, and b = 0. 



II. KAGOME TIGHT-BINDING SPECTRUM IN ZERO FIELD 

The kagome lattice is a two dimensional periodic array of corner-sharing triangles with 
three sites per unit cell. As illustrated in Fig. 1(a), a is the triangle edge length, A,B and 
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C denote the three sites in a unit cell . The unit cell can be taken as a rhombus of side 2a 
T and f 



with angles ^ and ^ at its vertices. In the reciprocal space k = {k x , k y ), the first Brillouin 



zone is a hexagon with a side length of Some points on the edge of the first Brillouin 
zone are shown in Fig.l(b): F(^,0) and"^^, and H(0, g). 

We first consider the tight-binding model for a particle making nearest-neighbor hops 
on the kagome lattice in zero magnetic field. The hopping matrix element is t > 0. This 
problem is readily diagonalized. Since there are three sites per unit cell, there are three bands 
of simultaneous eigenstates of the tight-binding Hamiltonian and the crystal momentum: 

ei(k) = e - 2t, 

,1/2 (2.1) 
e 2j 3(k) = e + t ±t 1 + 4 cos 2 (k y a) + 4 cos(k y a) cos(v3k x a) , 

where eo is the on-site energy of the orbitals. 

For simplicity, let us set t — 1 and e = 0. The spectrum of the second two bands (e 2 , e 3 ) 
ranges over the interval [—2, +4] and the edges of this interval are reached for k = 0. Close 
to k = 0, these two bands are parabolic with an isotropic effective mass. This effective mass 
is negative m* = —h 2 /(2ta 2 ) at the top of the upper band at e = 4. These two dispersive 
bands also touch at "Dirac points" at all the corners of the first Brillouin zone (such as G 
and H in Fig. 1(b)), at energy e 2 = e 3 = 1. 

The first band is completely non-dispersive at energy e 1 = —2. Thus the ground state 
of this particle is highly degenerate. Note also that the bottom of the lower dispersive band 
is also degenerate with this flat band at k = 0. Thus here at k = and e = —2 we have 
two degenerate bands of precisely the type discussed in the previous section of this paper. 



Since one band is flat, we have do = y^i + d 2 , m the notation of the previous section. The 
dispersive band, to quadratic order in the momentum, has energy e(k) — e(0) = tk 2 a 2 , so 
we have d = ^d 2 + d\ = ta 2 /(2h 2 ). 

Some of the degenerate ground states with e = — 2 in zero magnetic field are shown in 
Fig. 2. Fig.2(a) and (b) show the states at the Brillouin zone corners, while Fig.2(c) and (d) 
show the two degenerate states at the zone center with k = 0. For small magnetic field, the 
lowest Landau levels are composed of linear combinations of smooth envelopes times these 
latter two k = states. The two states shown are the eigenstates of the band operator <t 3 . 
By looking at their wavefunctions within a unit cell (one triangle) we can see that 03 does 
indeed measure an "internal" angular momentum that has the two possible values raj = ±1 
in this case, as anticipated above. 

We can further anticipate the results of applying a magnetic field by considering just 
these states. In the absence of coupling between the bands (i.e., for d± — 0), the states at 
the bottom of the bands are Landau quantized in the standard fashion, and their energy 
thus increases linearly with \B\. On the other hand, the degeneracy between the bands 
is split by the matrix element of the perturbation (di) that couples them. In the present 
case this perturbation is also linearly proportional to \B\. Thus whether the lowest-energy 
states increase or decrease in energy depends on the relative magnitudes of the coupling 
matrix element and the effective band masses. In terms of currents, a charged particle 
(here, a Cooper pair) in a given one-band Landau level has a quantized diamagnetic dipole 
moment that is independent of field. Its energy in a field thus increases with \B\. For this 
specific kagome lattice case, the wavefunctions within one unit cell also have currents and 
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dipole moments as in figure 2 c and d. When these degenerate states are coupled by the 
perturbation they are admixed and this can produce paramagnetic currents. The sum of 
this and the Landau diamagnetic current can be either diamagnetic or paramagnetic for the 
lowest-energy state depending on the relative strength of the parameters. In this case we 
find net paramagnetism in the ground states. 



III. SPECTRUM IN A MAGNETIC FIELD 

In the presence of a magnetic field B, the hopping term of the tight-binding Hamiltonian 
is modified by phase factors from the vector potential A, 20 , i.e., t — > te l7iJ , where 7^ is the 
phase factor between sites i and j: 

7 « = ^ / A • dl , (3.1) 

where $ — hc/e* is the flux quantum. 

The (tight-binding) Schrodinger equation in a magnetic field is thus the finite-difference 
equation: 

e<^ = i£e^, (3.2) 
3 

where the sum is over all the nearest neighbors of site %. It is convenient to measure the 
magnetic field in units of the flux quantum per elementary triangular plaquette of the kagome 
lattice. The flux through one triangle is $ = ^Ba 2 . Thus the "filling ratio", namely the 
fraction of a flux quantum through each triangle is 

$ _ y/3Be*a 2 
; "$o" 4 he ■ [6 ' 6) 

The flux through an elementary hexagon is six times this, 6/$ 0) while the total flux through 
a unit cell of the kagome lattice is 8/$ . 

When / is a rational number that can be written as 8/ = k/n, where k and n are 
integers with no common factors, the total flux through n unit cells of the kagome lattice is 
an integer multiple of $o- Then with this magnetic field, a gauge can be chosen so that the 
Hamiltonian has a discrete translational invariance with a unit cell containing n unit cells 
of the kagome lattice, and thus 3n sites. The spectrum of the Hamiltonian then consists of 
3n bands. These can be obtained numerically for any such rational / (for n not too large); 
an implementation of this is detailed in Yi Xiao's Ph.D. thesis. 21 The resulting spectrum, 
namely the kagome version of "Hofstadter's butterfly", is shown in Fig.3. 

The spectrum has various symmetries: It is only the flux modulo the flux quantum 
that affects the spectrum, so if / is changed to f + j with j any integer, the spectrum 
is unchanged. The spectrum is also unchanged on changing / to — /, because if if) is an 
eigenstate with energy e for field /, then its complex conjugate if)* is an eigenstate with the 
same energy e for field — /. These two symmetries are not special to the kagome lattice. 
The third symmetry is that if / is changed to / + |, this is the same as changing t to —t, so 
the spectrum is inverted. Thus the highest energy states for field near f — \ are equivalent 
to the lowest energy states near zero field. We will look at the behavior near zero field, but 
because of this last symmetry, what we find also applies near f — \- 
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IV. LANDAU LEVELS IN THE KAGOME LATTICE NEAR ZERO FIELD 



For small magnetic field, the spectrum of the kagome lattice tight-binding Hamiltonian 
shows standard equally-spaced Landau levels near the top of the bands, as shown in Fig.4. 
Here there is just the one band in zero field, with negative effective mass m* = —h 2 /(2ta 2 ). 
The resulting Landau levels, to linear order in the field, should therefore have energy 

e = t(A- 1 -^\f\(n+ 1 -)). (4.1) 

The numerical results for this portion of the spectrum are shown in Fig.4 for / > 
The first few Landau levels are clearly seen, and the asymptotic slopes at small / given by 
Eq.(4.1) are shown for comparison for the first 5 Landau levels. At these values of / the fit 
is not wonderful, but it does seem to be improving with decreasing /, as it should. 

More interesting is the bottom of the zero-field spectrum, where there are two degenerate 
bands. This results in non-standard Landau levels when the field is turned on, as discussed 
above. An expanded view of this part of the spectrum is shown in Fig.5. It shows some 
Landau levels whose energies are increasing linearly with field, as is normal, but these levels 
are not equally-spaced. In addition there are many states whose energies decrease as the 
field is increased. These latter states are the unusual feature we focus on here. Note that 
because of the symmetry of the spectrum discussed above, the same features appear at the 
top of the spectrum near the frustrating value of the magnetic field of / = 1/2. 

In the one-band case, knowing the effective mass in zero field and the charge is enough to 
determine the Landau level spectrum to linear order in the magnetic field. For the two-band 
case, on the other hand, (eq. 1.11) there is possibly the Bh ■ a term, which vanishes at 
zero field, so to determine the values of the three parameters making up b for the kagome 
tight-binding model, we cannot rely only on our knowledge of the zero-field spectrum. So 
we will now treat the low-momentum behavior to lowest order in the field. What we will 
find is that b = and d 3 = 0. 

The unit cell of the kagome lattice contains 3 sites around a triangle. We choose a set of 
3 basis states within this unit cell that respect as much as possible the three-fold rotational 
symmetry of this system about the center of the triangle. The wavefunction of a state 
is a complex number defined at each lattice point, but only its gauge-invariant properties 
are physical. These are the magnitude of the wavefunction, and the gauge-invariant phase 
differences between adjacent lattice points. For any wavefunction defined at the corners of 
the triangle in one unit cell, the total guage-invariant phase difference added up stepping 
the three steps around the triangle is 27r(/ + mj), where rrii is an integer. Note m 8 measures 
the state's "internal" angular momentum. The gauge- invariant phase difference between 
two adjacent lattice points is only well-defined modulo 2n, so let us restrict it to lie in the 
interval [— 7r, 7r]. Thus the three of these phase differences added up must be in the range 
[— 37r, 37r], and for the small fields / we consider here, this means the total guage-invariant 
phase difference around the triangle has to be near one of the three values: — 2tt, 0, or 2tt 
and mi is either -1, or 1. The basis states we will use are those where the guage-invariant 
phase differences are the same along all edges of the triangle, being 2n(f + mj)/3, while the 
amplitudes are identical at all three points. In the limit of small momentum and zero field, 
the upper band near energy 4t is made out of the m 8 = states while the lower bands near 
energy — 2t are made out of the two = ±1 states. 
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We take the tight-binding Hamiltonian and rewrite it in terms of these new basis states. 
This step is done in a specific Landau gauge. We then take the low-field, low-momentum 
limit, assuming that the (complex scalar) amplitudes of each of the basis states vary little 
between adjacent unit cells, so a gradient expansion can be made. The resulting three-band 
Hamiltonian can be written in gauge-invariant form as the matrix 

-2 + n 2 + f -iVS(U x -iU y ) (U x -tU y ) 2 
H = t\ iV3{U x + ill y ) 4 - 2IT 2 iV3{U x - iU y ) I . (4.2) 



IT, + illy) 2 -iV3(U x + illy) -2 + n 2 - 



B 



in units where a — h — 1 and <J>o — 2vr and where the rows/columns refer to = +1, 0, — 1 
in that order. Here we have gone to the orders in n and B that contribute to the spectrum 
to linear order in \B\ and quadratic order in n. Note that this continuum three-band Hamil- 
tonian conserves total angular momentum. To reduce this to the two-band Hamiltonian for 
the lowest bands only, we treat the coupling to the upper band in second-order perturbation 
theory. This eliminates the terms linear in B in the diagonal operators in the above matrix, 
and the two band Hamiltonian is precisely of the form expected (1.11) with b = 0, d 3 = 
and do = d\ = ta 2 /(2h 2 ). Finally, the Hamiltonian to this order is simply 

H = t[-2 + ^l[a (2a^a + 1) + 2(a_(at) 2 + a + a 2 )] (4.3) 
v3 

for / > 0. This Hamiltonian conserves total angular momentum m = m + mi. The resulting 
Landau level spectrum to linear order in the magnetic field is 

e± = t(-2 + ^=|/|(m + \ ± Vm 2 + m + l)) . (4.4) 

For m — — 1 and m = there is only the one state with m 8 = —1, because m cannot be 
negative. (There are no states with m < —1.) These states' energies are given by the plus 
sign in the above equation; they are the two lowest standard one-band Landau levels made 
out of only the = — 1 band. Since they are standard one-band Landau levels, they show 
diamagnetic currents. For m > 0, on the other hand, both bands enter and the eigenstates 
are linear combinations of the (m+l)th Landau level in the = —1 band and the (m — l)th 
Landau level in the = 1 band. The interaction between the two bands shifts the energies 
strongly, with the lower energy state being pushed to below energy —2t. The lowest energy 
state of all is the lower energy (minus sign in the above equation for e±) state with m — 1, 
which is at energy 

e mm = -t(2 + 8tt|/|(1 - ^)) - -t(2 + 3.37|/|) . (4.5) 

This state is a linear combination of the lowest Landau level (m = 0) in the = 1 band 
and the third Landau level (m = 2) in the = —1 band. And it is the interference 
between these two bands that allows this state to be paramagnetic, so its energy decreases 
with increasing \B\. 

The straight lines in Fig. 5 correspond to the lowest few of the upper (+ in Eq. 4.4) 
set of Landau levels. As / is decreased, the numerically obtained bands fit well to these 
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expected Landau levels. A more detailed plot of the lower bands below e = — 2t is shown 
in Fig. 6, in which this range of the spectrum is calculated for filling ratio from / = 1/16 
down to / = 1/256. For a rational filling ratio / = l/8n, where n is an integer, these lower 
bands consist of n — 1 subbands. Counting the states, we find that in the limit of small 
magnetic field 1/3 of all the states are in this set of bands below e = — 2t, corresponding to 
the number of states in the flat e = — 2t band at zero field. The lowest few Landau levels 
are also indicated by the straight lines in Fig.6, and again the numerically obtained bands 
fit well to these as / is reduced. 

To examine the very lowest Landau level in more detail, we carried out further numerical 
calculation of the bottom edge of the spectrum to / as low as 1/840. These results are shown 
in Fig.7 on log-log scales. The straight line corresponds to the small / behavior, Ae = 3.37/, 
derived above. As can be seen, it gives an excellent fit to the numerical results in this low 
field limit. 

This behavior in the kagome lattice is certainly different from the low field Landau levels 
of most other periodic lattices 1,5-7 , where the ground state energy increases as the magnetic 
field is turned on. In the case of the kagome latice, the presence of two degenerate bands in 
zero field allows these unusual Landau levels to occur. 



V. APPLICATION TO SUPERCONDUCTING WIRE NETWORKS 

The network equation for the mean-field phase boundary of a superconducting wire 
network, derived from linearized Ginzburg-Landau equation, has the same form as that of 
the tight-binding Schrodinger equation in the same geometry 4 . For a periodic network with 
equal bond length a, the superconducting network equation be can expressed as 

z cos j ■ ipi = ^2 ^je~ llij (5.1) 
^ i 

where ipi is the order parameter at site i, £ is the superconducting coherence length, z is 
the coordination number of the lattice, and the sum is over all nearest-neighbor sites of %. 
Compared to the tight binding equation of Eq.(3.2), the correspondence between the two 
system is given by 

e a . . 

- -<=>- 2cos - . 5.2 
t i 

The superconducting transition temperature of the network in a magnetic field is given by 
the upper edge (for t > 0) of the tight-binding spectrum. For the kagome lattice, z — 4. 
The upper edge of the kagome tight-binding spectrum is plotted in Fig. 8. 

For the phase boundary at the lower field limit (small /), our previous discussion of 
the low field Landau levels near e = 4 implies a linear relationship between the transition 
temperature and the magnetic field 

AT C (/) ~(|) 2 ~|/|, (5.3) 

where AT C (/) = T c (0) — T c (f) is the suppression of the critical temperature from its zero 
field value. 
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As shown in Fig. 8, strong maxima or minima in the superconducting phase boundary 
appear at magnetic fields / = p/8 for each integer p, which correspond to p magnetic flux 
quanta through each unit cell. At each of these fields, the tight-binding model is equivalent 
to a model in zero magnetic field with the hopping matrix elements within each unit cell 
being complex. Since there are three sites per unit cell, in each case there are three bands in 
the spectrum. When the magnetic field is changed slightly to / = | + A/, the small added 
magnetic field (of either sign) splits these bands into Landau levels, just like near / = 0. 
For all of these fields, except the special case / = 1/2, the upper edge of the spectrum is 
a state at crystal momentum zero in a nondegenerate band. Thus changes in the magnetic 
field produce conventional Landau levels like near / = 0. This causes the upper edge of the 
spectrum to decrease linearly in |A/|, resulting in a local maximum in the superconductor's 
transition temperature vs. field at those fields. The appearance of a local minimum in T c at 
f — 1/2 is due to the special degeneracy in the highest band that occurs due to the magnetic 
frustration at that field, resulting in Landau levels whose energy moves linearly away from 
the band center as the field is changed from this special value, as discussed above. 

A local minimum in T c vs. field also occurs at / = 1/2 in the so-called % or dice lattice 7 ' 8 , 
but for quite different reasons from the kagome lattice case. The dice lattice is the dual to 
the kagome lattice. Its elementary plaquettes are all rhombuses, and / is the number of 
flux quanta passing through each such rhombic plaquette. Like the kagome lattice, the dice 
lattice has three sites per unit cell. At / = 1/2 in this dice lattice, the spectrum consists of 
three dispersionless (infinite mass) bands 7 , with no degeneracy between the bands. For the 
dice lattice, the linear dependence of the upper edge of the spectrum on field near / = 1/2 
arises not from the formation of Landau levels, but from an "orbital Zeeman effect" . States 
with nonzero net magnetic moments of either sign and localized to just a few lattice sites 
can be formed out of the highest energy states at / = 1/2. The magnetic moments are due 
to currents circulating around the plaquettes and cause these states to have energies that 
depend linearly on the magnetic field. Although the kagome latice spectrum also has such 
a dispersionless band at / = 1/2, the localized states that can be made from this band in 
the kagome case have no net magnetic moment and thus no such orbital Zeeman effect. 
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Marc Schreiber, Yeong-Lieh Lin and Franco Nori for helpful discussions. This work was 
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FIGURES 




FIG. 1. (a) The kagome lattice is made up of triangular and hexagonal plaquettes with bond 
length a. Each unit cell has three sites, denoted here as A(x,y), B(x + ^a, y + \a) and 
C(x + ^a,y — \a). (b) The first Brillouin zone in reciprocal space. Some points on the edge 
of the first Brillouin zone: 0), G(^, £), and H{0, 
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FIG. 2. Some of the degenerate ground states with energy e = — 2 in zero magnetic field. 
00 = exp(i^) and uJ is its complex conjugate. The numbers next to the lattice points indicate the 
value of the wavefunction. The arrows along the bonds indicate the currents (supercurrents in the 
application to superconducting wire networks) and the signs at the centers of the triangles indicate 
the sign of the circulation of the currents. The k = states are (c) and (d). 
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FIG. 3. Spectrum of the kagome lattice tight-binding model in a uniform magnetic field. Here 
and in the following figures the energy is in units of t. The filling ratio / is the fraction of a 
magnetic flux quantum passing through each elementary triangle of the kagome lattice. 
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FIG. 4. Energy spectrum for Landau levels near e = At. The small / behavior of the levels 
corresponding to n = 0, 1, 2, 3, 4 are indicated here as solid lines. 
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FIG. 5. Lower subband structure close to e = —It. The calculated small / behavior of some of 
the Landau levels are indicated by the straight lines here. 
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FIG. 6. Subband structure below e = —2t for filling ratios down to / = 1/256. The calculated 
small / behavior of the 4 lowest Landau levels are indicated by the straight lines. 



17 



0.1 -. 



< 



0.01 - 



1E-3 




1E-3 



0.01 



f 



FIG. 7. Log-Log plot of the bottom edge of the spectrum close to e = — 2t, going down to 
/ = 1/840, with Ae = — 2 — (e m j n /i). The dots are from the numerics and the straight line 
indicates the expected small / behavior. 
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FIG. 8. Upper edge of the kagome lattice tight-binding spectrum for t > 0. The mean-field 
transition temperature T c of a kagome wire-grid superconductor is linearly related to e, so T c has 
a minimum at filling ratio / = 1/2. 
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